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Abstract 

Building on earlier results on holographic bulk dynamics in confining gauge 
theories, we compute the spin-0 and spin-2 spectra of gauge theories dual to the 
non-singular Maldacena-Nunez and Klebanov-Strassler supergravity backgrounds. 
We construct and apply a numerical recipe for computing mass spectra from cer- 
tain determinants. In the Klebanov-Strassler case, states containing the glueball 
and gluinoball obey "quadratic confinement", i.e. their mass-squareds depend on 
consecutive number as m 2 ~ n 2 for large n, with a universal proportionality con- 
stant. The hardwall approximation appears to work poorly when compared to the 
unique spectra we find in the full theory with a smooth cap-off in the infrared. 



1 Introduction 



One of the first successes of gauge/string duality was the computation of mass spectra of 
strongly coupled gauge theories from dual geometries [JJ [21 E] • These original papers were 
all concerned with A/" = black hole solutions, where it is hard to check the validity of 
the correspondence. Other work uses the Af = 4 superconformal theory (AdS bulk) plus 
a hard IR cut-off [I] that serves to imitate interesting field theory IR effects, but there is 
no running coupling. Studying Af = 1 super Yang-Mills theory coupled to matter seems a 
good compromise in many respects: these theories exhibit confinement, chiral symmetry 
breaking, running couplings and a rich set of mass spectra. 

For this "Non-AdS/Non-CFT correspondence" where the bulk is not AdS and the 
boundary field theory is not a CFT, much less is known about mass spectra, let alone 
correlators. A few brave attempts exist in the literature [HI El CD] , but as we explained 
in [TT], their results on mass spectra are at best inconclusive. In the present paper, we 
will compute the mass spectrum of the Af = 1 Klebanov-Strassler theory [12], a non- 
conformal deformation of the Af = 1 Klebanov-Witten theory [13] involving gluons and 
gluinos coupled to two sets of chiral superfields with a specific quartic superpotentialE] We 
also compute the mass spectrum of the Maldacena-Nunez background [16] . as a warmup. 
In this way, we will be able to address physical questions about how confinement works 
in these models, such as whether the theory displays "linear confinement" [32] . 

Despite this promising state of affairs, many authors have emphasized that the afore- 
mentioned theories are still quite far from real-world QCD. For instance, there are no 
open strings corresponding to dynamical quarks, hence no real meson spectra. On the 
good side, in [TBI EE], probe D-branes were added to the Klebanov-Strassler background, 
something which could be further developed using our methods. Real QCD is also non- 
supersymmetric, and we make heavy use of the existence of a "superpotential" W. How- 
ever, as indicated by the quotation marks, this superpotential is "fake" (in the sense of 
[2TJ]). but there are nonsupersymmetric examples where such structure exists irrespective 
of supersymmetry per se, see e.g. [21J. Last but not least, real QCD is not at large 't 
Hooft coupling A; this is very difficult to overcome with the present state of the art, but 
the recent progress in [22] at least shows that 1 / v^A corrections may not be unrealistic to 
obtain for these theories. 

Here, we take a step back from the effort to describe real QCD, and as mentioned 
above, focus on an example that is well-defined from a holographic point of view and 
see what can be understood about that theory. We introduce some new techniques, but 

■"■As far as we know, this expression was first used by Strassler at the Strings 2000 conference. It also 
featured prominently in the titles of Aharony's 2002 |5J and Zaffaroni's 2005 [3] lectures. 

2 More precisely, we consider only mass eigenstates that are dual to bulk solutions in the Klebanov- 
Strassler background, i.e. we do not consider meson mass spectra, which would arise from introducing 
flavor branes. Moreover, there is always an issue with contamination by Kaluza-Klein states in these 
theories, which we will ignore in the following. For nice introductions to the Klebanov-Strassler theory, 
see PUHi]. 
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Figure 1: Warp factors h(r). The string tension goes like /i -1 / 2 . From left to right 
the figures show the form of the warp factor in AdS (tension goes to zero in infrared) , 
Klebanov-Tseytlin (tension diverges at r s ) and Klebanov-Strassler (tension goes to 
finite value). 



conceptually the biggest difference from most of the literature is in things we do not do. 
We do not employ the "hard-wall" approximation (which amounts to taking AdS as in 
fig. [T]but with a finite IR cutoff). One argument that has been put forward to motivate 
the use of the hard-wall AdS model is that some physics should be insensitive to IR 
details. However, for computing mass spectra, a boundary condition must be imposed 
in the IR, so this question cannot really be insensitive to IR details (see e.g. [22] ). We 
discuss to what extent it is in section 15.71 Another common approximation, for example 
in the context of inflation [23], is the singular Klebanov-Tseytlin background [24] , as in 
fig. [H We insist on using the full Klebanov-Strassler solution and imposing consistent 
boundary conditions, and find that the commonly used approximations would likely not 
have led to correct results even for this theory, let alone for QCD. It seems reasonable to 
ask for explicit computational strategies to be developed concurrently with the search for 
real QCD, and we believe we have made progress on such strategies here. 

It is important to recall that the Klebanov-Strassler theory does not have a Wilsonian 
UV fixed point, but one can impose a cutoff in the UV to define the theory. One can think 
of this in at least three ways: 1) This is an intermediate approximation and the theory 
will be embedded in a more complicated theory with a UV fixed point like in [2H]; 2) we 
will glue the noncompact dual geometry onto a compact space that naturally provides a 
UV cutoff like in [26]; or 3) we will want to match to (at least lattice) data at that point, 
and beyond that the theory is in any case not asymptotically free (for recent ideas about 
this, see e.g. [27]). As far as this paper is concerned, any of these points of view may be 
adopted. 

Now to the new results in this paper. We expand upon our earlier work [TT] and 
present a numerical strategy, the "determinant method", to calculate the spectrum of 
regular and asymptotically subdominant bulk fluctuations. This puts us in the interesting 
position of being able to compute the mass spectrum unambiguously, but not being able 
to say what the composition of each mass state is without further information^ This 

3 As will become clear, this is also true to some extent in standard AdS/CFT, but often ignored. 
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information should ideally come from data, i.e. one should compute mixings at a specified 
energy scale where one has some control, and use the dual theory to evolve into the deep 
nonperturbative regime. To illustrate this, let us take the simpler example of a theory 
that does have a UV fixed point. The glueball operator O g , with conformal dimension 
A = 4, and the gluinoball operator Og, with A = 3, have a diagonal mass matrix in the 
conformal theory. As soon as we let the theory flow away from conformality, there is a 
nonzero correlator (O g Og). We can still diagonalize at any given energy, and the mixing 
matrix will be energy dependent. For theories that are not conformal even in the UV, 
there is no preferred basis labelled by the A eigenvalues. One can still contemplate an 
approximate labelling by A [in the KS theory, this corresponds to expanding in the ratio 
P = (number of fractional D3-branes)/(number of regular D3-branes)], but it is not yet 
completely clear how to implement this in the dynamics; as we pointed out in [11], the 
limit P — > does not commute with the UV limit one needs for the asymptotics. We will 
give ideas about this in sections 12.31 and 15.61 but we will not resolve it completely. 

The spectra we find have some simple features. The states of the Klebanov-Strassler 
theory come in towers, each of which shows quadratic confinement for large excitation 
number (i.e. m? ~ n 2 , where n denotes the excitation number within the tower). This 
confirms the claims of [T7j also for the spin-0 states of the KS theory. However, for low 
excitation number (i.e. roughly for the first two excitations in each tower) the structure 
of the mass values is richer. 

For the Maldacena-Nunez background, the spectrum we find is rather different. It has 
an upper bound, in agreement with the analytical spectra that we reported in [IT]. 

Let us also mention that there are other approaches to holography in this type of 
models, such as the Kaluza-Klein holography of [2H] ■ Some issues may appear in a different 
light in that framework, and that would be interesting to know. 

The organization of the paper is as follows. In chapter [2] we outline the general 
theory of bulk fluctuations as well as the correspondence between their spectrum and the 
spectrum of mass states in the dual gauge theory. The presentation is done such that it is 
applicable also in the case of non-asymptotically AdS bulk spaces. Based on this general 
material, a numerical strategy to calculate the specrum will be developed in chapter [3l 
As a warmup, our first application will be the Maldacena-Nunez theory [16] in chapter HJ 
before we come to the analysis of the mass spectrum in the Klebanov-Strassler cascading 
gauge theory in chapter [5J Some conclusions can be found in chapter El and many of the 
more technical details have been deferred to the appendices. 

2 Holographic mass spectra 

We want to calculate mass spectra of confining gauge theories from the dynamics of fluc- 
tuations about their supergravity duals. In this section, we review and develop the main 
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theoretical tools that are necessary for such a calculation. We will start, in Sec. 12.11 with 
a review of the treatment of the dynamics of supergravity fluctuations developed by us 
in [TTJ , generalizing a similar approach to the bulk dynamics in AdS / CFT [2UI I3TH |3"T] . 
Then, in Sec. \2.2\ we identify the bulk duals of massive states in the gauge theory as 
eigenf unctions of a second-order differential operator satisfying suitable boundary condi- 
tions. Most of this material belongs to AdS/CFT folklore, but we will fill in some details 
not commonly found in the literature. Finally, in Sec. 12.31 we will attempt to justify the 
identification of the gauge theory mass spectrum with the spectrum of bulk eigenfunc- 
tions by studying the pole structure of holographic two-point functions. We hope that this 
goes some way towards a formulation of holographic renormalization that holds beyond 
asymptotically AdS bulk spaces. 



2.1 Linearized bulk dynamics 

Let us start by reviewing the methodology we employ for studying fluctuations about bulk 
backgrounds that are dual to confining gauge theories. We use consistent truncations of 
Type IIB supergravity that describe a subsector of the 10-dimensional dynamics in terms 
of a 5-dimensional effective theory. Details of this trunction can be found in our earlier 
paper [UJ. Throughout this paper, we stick to a 5-dimensional effective bulk theory, but 
generalization to arbitrary dimension is straightforward. The effective theory is a non- 
linear sigma model of scalars coupled to 5-dimensional gravity, with an action of "fake 
supergravity" type, 

1 „ 1 



s = J dVs 



--R+ -G ab (<l ) )d,<p a dy b + V{<t>) 



(2.1) 



The "fake supergravity" property ^ implies that the scalar potential V(4>) follows from a 
superpotential W(<f>) by 

V{4>) = \G ah W a W h - V 2 . (2.2) 

where we used the shorthand W a = dW/d(p a . 

Background solutions, also known as Poincare-sliced domain walls or holographic 
renormalization group flows, are of the form 

ds 2 = dr 2 + e 2A(r) THi dx l dx j , 

13 (2.3) 

where the background scalars, 0(r), and the warp function, A(r), are determined by the 
following first-order (BPS-like) differential equations, 

d r A(r) = ~WU) , 
_ 3 KYJ ' (2.4) 

d r 6 a (r) = W a = G ab W h . 



4 Here "fake" does not mean that the system is necessarily nonsupersymmetric — the systems we 
considere here are supersymmetric — just that the formalism is applicable more generally. The relation 
between supergravity and fake supergravity was explored in |32l 133] . 
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Fluctuations about the Poincare-sliced domain wall backgrounds including the scalar 
and metric fields are best described gauge- invariantly. We will quote results from 
The physical fields are vectors in field space and scalars in spacetime, denoted by a a , plus 
a traceless transverse tensor, e*, where i,j = 1,2,3,4 are the indices along 4-dimensional 
Poincare slices. These fields satisfy the following linearized field equations, 

W a W h 8 \ / h h W h W c \ 0A 

5 a b D r + W a ]b - - -W5 a b \ U b c D r - W\ + ) + 5?e- 2A □ 



and 

-2A 



= o , 

(2.5) 

d'f. - ~Wd r + e~ 2A □ I z) = . (2.6) 

In ( 12.51) and (12.61) . □ = rfWidj, which will be replaced by — k 2 when we work in 4- 
dimensional momentum space. Also, we use a sigma-model covariant notation, lowering 
and raising indices of sigma-model tensors with G a b and its inverse G ab , and a sigma-model 
covariant derivative D a that acts as 

DaVb = ^a\b = d b if a ~ Q C ah ^c , (2.7) 

where Q c ab is the Christoffel symbol derived from the metric G a t,. Furthermore, we define 
a "background-covariant" derivative D r (r is the radial coordinate of (12.31) ) acting as 

D r ip a = d T ip a + Q\ c W h ip c . (2.8) 

So, our task is to solve (I2.5P and (12. 6p for the backgrounds of interest, imposing 
appropriate boundary conditions that we now describe. 



2.2 Bulk duals of mass states 

In AdS/CFT, the usual definition of the bulk dual of a boundary gauge theory mass state 
is a solution of the (linearized) bulk field equations that is regular and integrable. In the 
following, we will state these notions more precisely and present them in a form that is 
suitable also for application in non-AdS/non-CFT. 

As prototype for the equations of motion (12.51) and (12.61) . we consider a system of n 
coupled second-order differential equations of the generic form 



D 2 r - ^WD r + M(r) - e~ 2A ^ k 2 



, (2.9) 



where D r is the background-covariant derivative (12.71) . M is a symmetric field-space tensor 
(M ab = M ba ), and field indices have been omitted for the sake of brevity. We also assume 
that the radial coordinate r is defined in a domain tir < r < oo, where r w tir corresponds 
to the deep interior (IR) region, and large r to the asymptotic (UV) region of the bulk@ 



3 Wc allow ?'ir = — oo, like for example in a pure AdS bulk. 
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Let us start with the regularity condition. As we want the 10-dimensional configuration 
obtained via uplifting of an effective 5- dimensional solution to be singularity- free, we must 
impose that also the 5-d solutions be regular in the bulk. For nonsingular backgrounds, 
(12. 9p typically has a regular singular point at and no others for finite r j^| Regularity in 
the bulk, therefore, means regularity at riR. To make sure that our regularity condition is 
invariant under changes of field space variables, we require that the norm of the fluctuation 
vector is finite at r = r T R, 

Regularity condition: (a a G ab a b ) r=rm < oo . (2-10) 

We emphasize that one cannot simply demand regularity of the components of a a at 
r = tir, because r = can coincide with a coordinate singularity in the sigma-model 
metric G a b, when evaluated on the background. In the examples we consider, this indeed 
happens. 

From the 2n independent solutions of ( 12.91) . the condition (12.101) typicallyEl selects 
precisely n independent regular solutions, which we will denot^] by a rcgi ;, and n singular 
solutions, denoted by a S i ng ^. A generic solution of (12.91) is a linear combination of these, 

a(r) = d a regii (r) + q a singii (r) (2.11) 

with constants q and Cj. Thus, ( I2.10p implies that 

5i = for all i. (2.12) 



Now, consider the integrability condition. Let us take a small detour and consider the 
bulk Green's function G(r,r'; k 2 ), which satisfies 



D 2 - ^WD r + M(r) - e~ 2A ^ k 2 



G(r, r'; k 2 ) = - e~ 4A ^ S(r - r') , 



(2.13) 



where the factor e 4A on the right hand side is the metric factor 1 / ^fg from the covariant 
delta functional. The Green's function can be written in terms of a basis of eigenf unctions, 



G(r,r',k 2 ) = J2 



a\{r)a x {r 
k 2 + m x 



2 ' 



(2.14) 



where the functions a.\ satisfy (I2.9P for k 2 = —m\. (Again, we omit the matrix indices, 
and the indices of the two cia's are not contracted.) Substituting ( 12. 14ft into (12 . 131) yields 
the completeness relation 



ax(r)a x (r') = e" 2A(r) 5(r - r') 



(2.15) 



A 



6 By "typically" we mean that this is the case in all regular holographic configurations known to us. 
It would be interesting to find more precise conditions for these statements. 

7 We hope there will be no confusion between the index i — 1, . . . , n and the field space index a. 
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from which one can deduce the orthogonality relation 

J dr e 2A ^ a x (r) ■ a a (r) = 6 Xa . (2.16) 

With the dot product we denote the covariant contraction of indices. Eq. (12.161) provides 
the condition for the eigenstates a\ to be integrable. Due to the factor e 2A (for general 
dimension d it would be e^ -2 ^), the integral measure is not the covariant bulk integral 
measure that one might have expected. 

Now, the spectrum of bulk eigenfunctions a a is identified with the mass spectrum of 
operators dual to the fields a. (We will also say that the eigenfunction a.\ is the bulk dual 
of a state of mass m\.) Clearly, this spectrum depends on the boundary conditions one 
imposes. One boundary condition is given by the regularity condition (12.101) . the other 
follows from (I2.16p . namely 

Integrability condition: J dr e 2A ^ o(r) • a(r) < oo . (2-17) 

For large r, i.e., in the UV, (12.91) has 2n independent asymptotic solutions. Let us 
denote the n dominant ones by cidom.i, and the n subdominant ones by a su b,i (by "domi- 
nant" we mean leading in r). Again, a generic solution of (12.91) can be written as a linear 
combination of these, 

a(r) = di a sub ,i(r) + d~i a domii (r) , (2.18) 

with constants di and di. 

It can be checked in the various cases we consider that the asymptotically dominant 
behaviours, ad om ,i, are not integrable in the sense of (12. 17ft . whereas the subdominant 
behaviours are integrable. Thus, (12.171) is equivalent to 

di = for all i. (2.19) 

To summarize, the eigenfunctions Oa, interpreted as the bulk duals of boundary field 
theory states with mass m^, are such that 

c i = d i = i = l,2, ...n, (2.20) 

from which follows that 

a\{r) = c X a a reg ,j(r) = d Xii a subii (r) . (2.21) 

2.3 Pole structure of holographic 2-point functions 

The masses of states (particles) in a quantum field theory manifest themselves as poles 
in the two-point functions of operators, if there is a non-zero probability that these states 
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are created by the operators from the vacuum 



d A x e ikx (O 1 (x)O 2 (0)) = J2 



(O|0 1 (O)|A)(A|0 2 (O)|O) 
k 2 + m\ 



+ c.t. . 



(2.22) 



For comparison with this general formula, we will now try to obtain the pole structure 
of holographic 2-point functions from the linearized dynamics of the bulk fields that we 
described in Sec. 12.11 But we must start with a disclaimer. As we have not systemati- 
cally dealt with the renormalization and dictionary problems (see the discussion in [IT]), 
our results will depend on two strong assumptions, motivated by how things work in 
AdS/CFT. A better understanding of these assumptions as well as the meaning of cases 
in which they are not satisfied would be very desirable. 

Let us start with the asymptotic expansion of the bulk fields f)2.18p . For a given bulk 
field a, the coefficients di and d{ clearly depend on the choice of basis functions a su b,i 
and cidom,i, respectively. In order to remove some of the ambiguities, we assume that the 
term containing k 2 in the equation of motion ( 12.91) is asymptotically suppressed, so that 
the leading terms of the asymptotic solutions are independent of k 2 , and the sub leading 
terms can depend only on powers of k 2 (assumption Then, we can choose a basis of 
asymptotic solutions such that each has a distinct leading behaviour. These solutions are 
classified according to their asymptotic growth with the radial coordinate into dominant 
and subdominant solutions. With each dominant asymptotic solution adom,«, we associate 



of basis functions, but we now argue that these correspond to field theory ambiguities. 
First, the normalization of the leading terms in the asymptotic basis solutions can be freely 
chosen. In field theory this corresponds to a choice of normalization of the operators Oi, 
which always drops out in physical scattering amplitudes. Second, we have the freedom 
to add to a given dominant solution multiples of asymptotic solutions of equal or weaker 
strength, with coefficients polynomial in k 2 . (Solutions of equal strength can be added 
only with coefficients independent of k 2 .) This is also expected from the dual field theory, 
since operators of a given dimension can mix with operators of equal or lower dimensions 
(higher dimension operators being suppressed by our large UV cutoff) which implies that 
generically, there is no unique way to define the renormalized operators. Therefore, the 
remaining ambiguities in the choice of an asymptotic basis of solutions reflect the usual 
freedom in the definition of field theory operators. 

For what follows, we do not need to make particular choices for the subdominant basis 
solutions, a su b,i- This will be clearer once we obtain our final result of this section, i.e. 



8 Typically, terms that are formally infinite and analytic in k 2 (contact terms, "c.t.") are needed to 
make the sum over the spectrum convergent. The finite parts of these countcrtcrms are renormalization 
scheme dependent. For a simple example, see Appendix |A"1 

9 This assumption is quite strong, and as we shall see, it holds in the KS system, but not in the MN 
system. 

10 This does not make reference to components of the field space vector a(r), which would require the 
operator O to carry a field index a and make it transform under bulk field redefinitions. 




Two ambiguities still remain in the choice 
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Now for the second assumption. In a field theory calculation of correlation functions 
of the operators O i: one adds a source term to the Lagrangian0 

AC = -diOi . (2.23) 

In holography, the sources are identified with the coefficients di of (12.181) . whose depen- 
dence on k 2 was suppressed in that equation. The second strong assumption we make 
is the generic form of the exact one-point function of the operators Oi. In AdS/CFT, 
the response functions di of (12.181) encode the exact one-point functions, up to scheme 
dependent terms that are local in the sources [3U EHJ [36l [37] o In view of the ambiguities 
discussed above, we assume 

(CM^Lact = ^j(k 2 ) dj{k 2 ) + local terms , (2.24) 

with a matrix Yy(/c 2 ), which we will be able to determine further below (assumption 
2). Part of the assumption is that the poles of the connected 2-point functions arise in 
dj, and that Yij does not give rise to additional poles. In AdS/CFT, (12.241) follows from 
holographic renormalization, but unfortunately we have no proof of it in this more general 
setting. 

We will now derive the general pole structure of holographic 2-point functions. To 
start, consider the general formula for a solution a(r, k 2 ) of (I2.9P in terms of the Green's 
function and prescribed boundary values. Let ro be a (large) cut-off parameter determin- 
ing the hypersurface where the boundary values are formally prescribed. Remembering 
that neither the Green's function nor its derivative vanish at the cut-off boundary, we 
havd^] 

a (r, k 2 ) = e 4A ^ [{D ro a(r , k 2 )) ■ G(r , r; k 2 ) - a(r , k 2 ) ■ D ro G(r , r; k 2 )} , (2.25) 

where a(r , k 2 ) and D ro a(r , k 2 ) are the prescribed values of the field and its first derivative 
at the cut-off boundary, respectively. Since r is an unphysical cut-off parameter, we 
must ensure that the bulk field a(r) remains unchanged when ro is varied. This is easily 
achieved, if, together with a change of the cut-off, r — > r + Sr , the boundary values are 
changed by 

5a(r , k 2 ) = (D ro a(r , k 2 )) 5r , 5D ro a(r , k 2 ) = (D 2 a(r , k 2 )) 5r , (2.26) 



11 The minus sign is a useful convention. Let (Oi) t be the exact 1-point function (i.e., the 1-point 
function in the presence of finite sources dj). Then, a connected (n + l)-point function involving Oi is 
given by the n-th derivative of (Ci) exact with respect to the sources dj. With a plus sign, one would have 
an additional factor ( — 1)". 

12 Imposing the regularity condition (|2.10p . only one set of coefficients is independent, which is taken 
to be di. 

13 This formula follows from (|2.13|) upon multiplication by e 4A ^ r ' a(r) from the left, taking the integral 
over r, integrating by parts and using the field equation (|2.9[) . The IR boundary does not contribute, 
because e 4A vanishes there. The reason for this is that r = should correspond to a single point of the 
bulk space, which is only guaranteed if e 4A vanishes there, c.f. 
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and the second derivative, D 2 o a(r , k 2 ), is determined by the equation of motion ( 12. 9ft . 
To assure (12.261) . we determine the formal boundary values at the cut-off, o(r , k 2 ) and 
D ro a(r , k 2 ), from the generic asymptotic behaviour (12.181) . with coefficients di and dj 
fixed. After inserting (T2~Tg|) and (l2~T4"j) into (l2~2"oD . we obtain 



a ( r ,A; 2 )=e 4A ^^^ ,x 



4A(r ) \^ 

k 2 + m\ 

A (2.27) 
x {di [(D ro a suh j(r , k 2 )) ■ a A (r ) - a sub ,z(r , k 2 ) ■ D ro a x (r )] + 

+ di [(D ro a dom j(r , A; 2 )) ■ a A (r ) - o domjZ (r , k 2 ) ■ D ro a x (r )] } . 

To continue, we observe that for very large r , the term on the second line of (I2.27p . 
containing only subdominant solutions, is much smaller than the term on the third line. 
Therefore, we drop it. Moreover, as we are interested only in the pole structure, we con- 
sider k 2 very close to — m\ and expand the numerator keeping only the leading term, i.e., 
we replace k 2 by — m\ in the numerator. Finally, we use the fact that the eigenfunctions 
are purely sub dominant I* 4 ] 

a\(r) = d x ,i a sub ,i(r, -m\) . (2.28) 



This yields 



a(r, k 2 ) = d, Zy(-ro') a sub ,(r, -m\) , (2.29) 

fc 2 «— m x K -\- iTl x 



with 

Z^k 2 ) = e 4A « [(D r a dom ,(r, k 2 )) • a subj (r, k 2 ) - a domjl (r, k 2 ) ■ D r a Buh>j (r, k 2 )] . (2.30) 

Notice that, by virtue of the equation of motion (12. 9p . Zy is independent of r. 

Thus, after reading off the response function d f from (12.291) . we obtain the poles of the 
connected 2-point function, using ( 12.241) and differentiating with respect to the source dj, 

(0^0 J (-k)} = J2Yui~ml)Z 3 A-rnl)^^ + c.t. . (2.31) 

By symmetry, we find that Yy equals Z^ up to normalization, 

Yij(-ml) = N x Z ij (-m\) . (2.32) 

Therefore, defining also 

Z\,i = Ziii(—m\) d x>i i 

= e 4A{r) [(D r a dom:i (r, -m\)) ■ a x (r) - a doilM (r, -m\) ■ D r a x (r)] , (2.33) 



14 At this point one may wonder where the dominant part of o comes from. It arises from the sum over 
the spectrum in (|2.27p . in particular from the UV contribution. For the simple case of AdS bulk, we show 
this in Appendix lAl However, it does not contribute to the poles. 
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our final result is 



(0 4 (*)0i(-AO> = £>A 



A 




f + c.t. . 



(2.34) 



We note that the Z\,i are independent of the choice of the subdominant basis solutions, 
because the normalization of the eigenf unctions a a is fixed by (12. 16ft . 

The determination of the normalization factors N\ is an open problem, whose solution 
is intimately related with finding a proof of (12.24)) . We conjecture that N\ = 1, as can be 
easily checked for asymptotic AdS bulk spaces (see Appendix [AT). 



In the previous section, we characterized the duals of mass states as regular (in the IR) and 
asymptotically subdominant (in the UV) solutions of the linearized bulk field equations. 
To calculate the spectrum, one solves a system of second-order differential equations with 
boundary conditions specified at two endpoints. The standard numerical strategy for 
this is "shooting", in which one solves the corresponding initial value problem with two 
boundary conditions (field value and derivative) specified at the initial point, and then 
varies the initial condition trying to find the desired boundary value at the final point. For 
coupled multifield systems, this is a laborious process. Fortunately, since we are interested 
only in the spectrum and not in the explicit dual bulk solutions, we can adopt a much 
simpler strategy. It involves computing a single function, the determinant of a matrix, 
that depends on the momentum k 2 and the vanishing of which signals the presence of a 
mass state. 

In Sec. 13.11 we shall present the general idea of our strategy in its simplest implemen- 
tation involving the numerical solution of the field equations as an initial value problem 
starting close to riR. This simple implementation presents some numerical challenges, 
which we describe, and which are overcome with the refined method described in Sec. 13.21 
Further details of numerical issues will be deferred to Appendix [El 

3.1 Mass Spectrum: Determinant Method 

We begin by imposing regularity in the IR. This means that we consider the equations of 
motion for r pa r^, find independent solutions of these, and impose q = in the language 
of Sec. 12.21 Starting with these initial conditions, one numerically evolves to larger r to 
obtain the n regular bulk solutions a reg ,j(r). For large r, the asymptotic behavior of each 
of these solution can be expanded in terms of the UV-dominant solutions Odomj? so t na t 



for some fc 2 -dependent but r-independent constants 7^. (The order of indices has been 
chosen for later convenience. As a reminder, both a regi j(r) and ad om ,j( r ) are n-component 



3 Finding Mass States 




(large r) , 
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vectors in field space, but we are suppressing the field-space index.) Thus, a general 
regular solution (cf. (12.111) ) will behave as 



areg(r) « CiljiddamA 7 ") (larger). (3.2) 

In order for a reg to qualify as the bulk dual of a boundary field theory mass eigenstate, 
the remaining condition in (I2.20p is dj = for the fluctuation to be integrable. Thus, 

di = ^ijCj = for all % -<=>- mass state . (3-3) 

Regarding this as a set of n equations for the coefficients Cj, it is necessary and sufficient 
for the existence of non-zero solutions that the determinant of the matrix 7^ be zero: 

det 7ij(/c 2 ) = -<=>- mass state with m 2 = —k 2 , (3.4) 

where we have restored the dependence of the matrix 7^ on the 4- momentum k 2 . 
To compute 7^, let us make the field-space component index explicit in eq. (13.11) : 

Keg^n » 7i<(a<taa(*)) i (3-5) 

where a reg (<2) and adom(^) are viewed as n x n matrices, with column index % and row 
index a. Inverting this equation, the coefficient matrix 7^ can be found by calculating^! 

lij ~ [(a dom y l (a Icg )].. (3.6) 

for some large value of r, where the subdominant parts of sufficiently suppressed. 

We would like to emphasize that formula (I3.6P renders the coefficients 7^ independent of 
the coordinates in field space — the component index a is contracted. Hence, also the 
mass condition ( 13. 4ft is independent of the choice of sigma-model variables. 

However, as announced above, the IR — > UV determinant method does not work 
straightforwardly for coupled multifield systems in practice. The method is only suffi- 
ciently insensitive to numerical error if all dominant solutions ad om ,i are comparable in 
magnitude, such that the integration error from the stronger ones does not swamp the 
significant digits of the weaker ones. For 1-component systems, this is obviously never a 
problem, and it turns out also not to be a problem for the 2-component system in the 
MN background (Section 14.3ft . In general, and in particular for the 7-component spin-0 
system in the KS background, it is better to tread more carefully and only evolve to a 
prescribed midpoint value, as we will describe in the next subsection. 

One could have imagined that this numerical difficulty would be alleviated by evolving 
the other way (UV — > IR). This is actually partially true, but the method in the next 
subsection is still more useful (see Appendix [E] for a few more details). Indeed, imposing 
the relevant initial conditions at a large value of r, one can calculate numerically the n 

15 This equation is intended in the sense of matrix multiplication, so there is no G a b involved. 
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UV-sub dominant solutions a su b,i- Decomposing their behavior at r ps tir, using (12.111) 
into 



the same argument as before affirms that an a su b dual to a mass eigenstate exists, if 
det^/y = 0. The coefficients 7^ can be determined by matching the components of a su b 
to the generic singular behaviour. 

The build-up and growth of integration error on the way from the asymptotic region to 
rpR is an issue here and must be resolved. This can be done by making sure the UV cutoff 
is not too large. Then, the integration error, which generally grows like the asymptotically 
weakest solution, remains small. 

3.2 Mass Spectrum: Midpoint Determinant Method 

The midpoint method is numerically the most stable and combines the two previous 
approaches. One calculates the asymptotic solutions analytically both in the UV and in 
the IR. Then, the n regular solutions a reg ^ are evolved from r r IR up to a mid-point 
r mid , and the n subdominant solutions a su b,i are evolved from large r down to r mid . Then, 
one tries to find linear combinations of the regular solutions and of the subdominant 
solutions, respectively, such that their value and first derivative at r mi d match: 



As before, the necessary and sufficient condition for the existence of a mass eigenstate 
is given by the determinant of a matrix (this time 2n x 2n) being zero. Put differently, 
rather than demanding that a given individual field match smoothly across the midpoint, 
we allow for the situation that only some linear combination of the basis fields matches 
smoothly, and this linear combination is determined by the coefficients c and d. 

If the solutions we are matching were approximate analytical solutions and the mid- 
point were the classical turning point, this midpoint matching would be precisely the 
WKB approximation. Of course, our solutions will not be approximate analytical, but 
numerical and exact (up to integration error), so it is no approximation in that sense. 

4 Maldacena- Nunez system 
4.1 General Relations 

The effective 5-d model describing the bulk dynamics of the MN system contains three 
scalar fields (g, a,p) and is characterized by the sigma-model metric [381 E] 




(3.7) 





G ah d^ a d^ h = d^g + e~ 29 d^a + 2Ad^p , 
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and the superpotentia 

W = -\ e 4p [(a 2 - l) 2 e~ 49 +2(a 2 + 1) e~ 29 +l] 1/2 . (4.2) 

Let us summarize what Poincare-sliced domain wall backgrounds this system admits. 
In what follows, we shall denote the background values of the scalar fields by g, a and 
p, and for their fluctuations we will use the gauge-invariant sigma-model vector a = 
(8g, 5a, 5p) T The background equations ( 12. 4p are most easily solved in terms of a new 
radial coordinate, p, defined by 

d p = 2e- 4p d r . (4.3) 
We shall need only the explicit solutions for a and g, which are given by 

2p 



sinh(2p + c) 



e 29 = 4p coth(2p + c) - (a 2 + 1) . (4.4) 



The integration constant c, which appeared for the first time in [39], can take values 
< c < oo. The regular MN solution corresponds to c = 0, while c > leads to singular 
bulk geometries. The solution (14.411 is defined for p > p, with p defined by e 29 ^ = 0, from 
which one obtains the relation 

coth(2p + c) = - 1 . (4.5) 
Zp 

Hence, p takes values between and 1/4, with p = for c = and p = 1/4 for c = oo. 
For c > 0, p = p is the location of the singularity. 

For the scalar p we will only need its relation to the warp function A (see ( 12.31) ). which 

is 

e -2A e -S P = C 2 ^ (4 6) 

where C 2 is an integration constant setting the 4-d scale. For later convenience, we set 
C 2 = 1/4. 

Let us turn to the fluctuation equations (12. 5ft and ( 12.61) in terms of the radial variable 
p. After using the relations for the background, one straightforwardly finds that ( 12.51) 
becomes 

[(d p -M)(d p -N)-k 2 ]a = 0, (4.7) 
where we dropped field indices, and the matrices M and N are given by 

A ( W a W h 
N a b = 2e~ 4p { d h W a - 



M\ = -N a b - 4e" 4p (G a bc W c - W6$) 



16 We have adjusted the overall factor of the superpotential of [38 to our conventions. 
17 The superscript T denotes the transpose, because o is a column vector. 
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Note that M and N are independent of p. As was observed in [TT], the fluctuations of p, 
which are described by the component a 3 , decouple from the other two. This can be seen 
from (14. 8 p as follows. From (14. Xp and (14 .2p we find the identities 

W 3 = 4W , W 3 = -W , cW a = AW a , (4.9) 

6 

from which follow iV a 3 = N\ = 0. Moreover, one easily checks that Q 3 bc = Q a 3c = 0, so 
that M% = Ae-* p W5l, and M\ = Ae~^W5l. Thus, a 3 satisfies 

[{d p - 4e" 4p W)d p - k 2 ] a 3 = . (4.10) 

This is just the equation for a massless scalar field in the domain wall background. 
Similarly, (12. 6p gives rise to 

[(d p - 4e~ 4p W)d p - k 2 ] t) = . (4.11) 

which is identical to 114. lOR . 



4.2 Analytic solutions in singular background 

We will briefly introduce the analytical solutions for the fluctuations in the singular back- 
ground case c = oo derived in [TT]. These analytical solutions give useful indications 
for glueball and gluinoball masses that can be compared with the numerical analysis in 
Sec. 14.31 Further details can be found in our earlier paper [TT] . 

Let us start with (14.71) . In the case c = oo, the matrices M and N are very simple 
and diagonal, 

M = diag ( -8p + 4 - -, -2 + -, -8p 

AP ~) , } , P P ' (4-12) 

N = diag — , - - 2,0 

V p p 

With these expressions, ( 14. 7p can be solved analytically. We shall illustrate the analysis 
considering the independent solutions for the second component, 

a 2 ~ e-^' 2 > I $ p ~ f °' 2 5 Z) ' (4.13) 

where $ is the confluent hypergeometric function of the first kind (or Kummer function 
M), and a and z are defined by 18 l 

a = \{l + k 2 y\ z =\{p-l) ■ ( 4 - 14 ) 



18 The combination p — 1/4 is simply p — p; see the discussion around 
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The two solutions in ( 14. 13ft have different behaviour at the singularity, i.e., for p ps p = 
1/4. Notice that both solutions satisfy the regularity condition (12.101) . so we do not have 
the standard means of excluding one of them. In other words, in the background with c = 
oo, the regularity condition (12.101) does not uniquely specify the IR boundary value that 
should be imposed on the eigenfunctionsl^l Let us tentatively discard the solution with 
the leading small- z behaviour of ( 14.131) . which will be justified by the numerical results 
in the next section. This leaves us with the first solution. The generic asymptotically 
dominant UV behaviour is absent in this solution, if a takes one of the values pTT] 

a = l(4n + l) , n= 1,2,3,... . (4.15) 
6 

As the second component is dual to the gluino bilinear, these states may be interpreted 
as gluinoballs. 



Applying the same argument to the solutions of the first component [TT] , we find states 

for 

a = -Un + 3), n=l,2,3,.... (4.16) 
6 

The spectra (14.151) and (14.161) can be combined into 

a = l(2n + 3) , n = 1,2,3,... , (4.17) 
6 

with odd and even n corresponding to (I4.15P and (14.161) . respectively. Using (I4.14p . we 
can rewrite (I4.17P as 

9 4n(n + 3) 9 , 

m 2 = . K | ' =1- - —rz , n = l,2,3,.... (4.18) 



(2ra + 3) 2 (2ra + 3) 



We would like to remark that, in [IT] , we had obtained a spectrum identical to ( 14. 18ft , 
with the addition of a massless gluinoball for n — 0. The above heuristic argument, that 
the solution with leading small- z behavior should be discarded, disposes of the massless 
gluinoball. This choice will be justified a posteriori by the numerical results. 

As we move from the background with c = oo to the regular background with c = 0, 
we expect these values to change. Moreover, in the regular background, no heuristic 
argument is needed, as the regularity condition will exclude half the solutions. 

From the solutions of the third component, no discrete mass values are found |llj . 



4.3 Numerical solutions in the regular background 

Let us now consider (14.71) in the regular background, i.e., for c = 0. The matrices M and 
iV are considerably more complicated than in the previous case, and we list their entries 
in Appendix [B] Thus, the full equations of motion will be solved numerically. 

19 One can show that the same behaviour is found for all c > by expanding the matrices M and N for 
p ~ p. The regular case c = is qualitatively different from these, because the limits c->0 and p — ► p 
do not always commute. To see this, consider, e.g., da(p)/ dp at p = p. 
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For our purposes, the simple approach of Sec. 13 . II is sufficient for the MN solution. (It 
will not be sufficient for the KS solution.) Hence, let us first consider p in the vicinity of 
p = 0, and find the regular solutions of (14. 7p . Expanding the matrices M and N about 
p = yields 



M 



N 



/-f + J J-g o \ 

p a p n 
% i-5£ o 

3 p 3 

o o -~-f 

9 3p 15 



+ o (p 3 ) 



-4p 

V 



4 + f o 
o oy 



+ ^ (p 3 ) • 



(4.19) 



(4.20) 



Recall from the discussion at the end of Sec. 14.11 that the zeros in (14.191) and f)4.20p are 
exact. 



Using (14.191) and (14.20p , it is straightforward to solve (14 .7ft in terms of series expansions 
about p = 0. The independent regular behaviours are 



ai 



'1 N 

o] +o( P 4 ) 

v0> 



(4.21) 



a 2 = p 2 \3\+0(p*) 

,0; 



03 



i + i*v 




(4.22) 



(4.23) 



The other three solutions go as 1/p in the leading term and will be discarded. 

The asymptotic region (p ^> 1) is well described by the singular background considered 
in Sec. 14.21 Hence, a convenient basis of generic dominant asymptotic behaviours can be 
found from the analytical solutions of [TT]. Combining them to a matrix, where each 
column is an independent solution, we get 



' P 



fldc 



-l/4-3af/2 






P 





l/4-3a/2 





(4.24) 



Putting into practice the method of Sec. 13.11 we have searched for zeros of det 7, where 
the matrix 7 is defined by ( 13.61) . with cidom given by (14.241) . and a reg being the matrix of 
independent numerical solutions at large p. In order to calculate them, it is convenient 
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Table 1: Comparison between numerical results (m 2 ) and analytical estimates for 
the mass spectrum in the MN background. 
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to rewrite ( 14. 7ft system of first-order equations 



In order to calculate the three independent regular solutions, we impose the three initial 
conditions (I4.2ip - fl4.23p at very small p. 

Our results are as follows. For the system involving the first two components of a, we 
obtain discrete mass values. The first few of these are listed in Tab. [Hand can be compared 
to the estimates from the singular background, (14.181) . To refine this comparison, in Fig. [2] 
we plot the function 

/ (777,77) = (l -m 2 y 1/2 - -77 , (4.26) 

o 

which leads to the mass values 

m2 = 1 - <dW ■ (4 - 27) 

If (I4.18P were exactly true, f(m, n) would be identical to 1. We find that the numerical 
values consistently come out smaller than those from the analytical approximation, and 
that they cluster into two distinct spectra. For even n, / ~ 0.95, whereas for odd 77, 
except for the first few values, / ~ 0.78. This appearance of two spectra is not too 
surprising, if we remember that, in f)4.18p . values for even and odd 77 originated from 
the first and second component of a, respectively. These two components couple in the 
regular background so that their spectra will combine into one. However, the numerical 
results indicate that there are two different kinds of bound states in this spectrum. 



20 Notice that both a and b contain 3 components, so that (|4.25[) is a system of 6 first order ODEs. It 
is convenient to consider the third component of a separately, as it fortuitously decouples from the other 
two. Thus, we have a 4-component and a 2-component system of first order ODEs. 
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Figure 2: The function (|4.26p for the numerical mass spectrum in the MN back- 
ground. 



From the massless scalar equation ( 14.101) . which describes both the third component 
of a and the spin-2 field e*, we have not found any discrete mass states, confirming the 
analytical result of Sec. 14.21 



5 Klebanov-Strassler system 
5.1 KS Background 

The effective 5-d model describing the bulk dynamics of the KS system contains seven 
scalar fields. We will use the Papadopoulos-Tseytlin [38] variables (x,p, y, $, b, hi, h 2 ). In 
the following, we shall briefly summarize the general relations for this system and the KS 
background solution. For more details, see our earlier paper [TT] and references therein. 

The sigma-model metric is 

G ab d^ a d fl (J> b = dpx&x + Qd^p + -dpy&y + -d^VQ + ^ e*~ 2x d^bd%+ 

+ \^~ 2x b~ 2y d^h 1 -h 2 )d^h 1 -h 2 )+e 2 yd,(h 1 + h 2 )d^h 1 + h 2 )] , (5.1) 

and the superpotential reads 

W = ~ (e~ 2p - 2x + e 4p cosh y) + - e ip ~ 2x (Q + 2Pbh 2 + 2Ph 1 ) . (5.2) 

Here, Q and P are constants related to the number of D3-branes and wrapped D5-branes, 
respectively. 

It is useful to introduce the KS radial variable r by 

d T = e~^d r . (5.3) 
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Table 2: Mass spectrum of spin-2 glueballs in the KS background (m 2 < 100), and 
comparison with Krasnitz' WKB results [30]. Krasnitz' values have been normalized 
such that the seventh masses (the highest he calculated) agree. 
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7 


17.1 




7 


17.06 


14 


60.41 







Here and henceforth, a field variable denotes the background of that field (p in this case) , 

while the sigma-model covariant a will be used for fluctuations. In terms of r, the KS 
background solution of (12.41) is given by 

$ = $ , (5.4) 

e y = tanh(r/2) , (5.5) 

sinhr ' ^ 

hi = -^p + Pe*° coth r(r coth r - 1) , (5.7) 

^ r coth t — 1 . 

h 2 = Pe*° — , 5.8 

sinn t 

2 e 6 P +2x = CQth r r ^ ( 5 ^ 

3 sinh t 

e 2,/3-4 P = 2p2 e$0 3-2/3^) sinh 4/3 r ? (510) 



with 



OO 



/?9 coth 79—1 
dtf u [2 sinh(2tf) - 4tf] 1/3 . (5.11) 

sinh xJ 

T 

Moreover, one can show that the warp function A satisfies 

e -™-S P „ ( e -6 P -2, sinh r j2/3 (512) 

where the proportionality factor depends on an integration constant that sets the 4-d 
momentum scale. 



5.2 KS spin-2 spectrum 

Let us start the numerical analysis of fluctuations in the KS background with the simplest 
case: the free massless scalar describing the transverse traceless components of fluctuations 
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FIGURE 3: Mass-squared values of spin-2 states. The solid line represents the fit 

dsn. 

of the bulk metric. Specifically, we are interested in the spectrum of spin-2 glueballs, 
expanding on earlier work by Krasnitz [7J HQ] . 

The equation of motion (12. 6p . after substituting the appropriate relations from the 
previous subsection, takes the form 

(d 2 T + 2 e~ 6p - 2x d T - k 2 e~ 2A - 8p ) e = , (5.13) 

where we omitted tensor indices on c. Defining 

h(r) 



It) 



h(0) 



(5.14) 



and appropriately choosing the integration constant in the warp function (15.121) such that 



e -2A-8 P= ( e -6 P -2x sinhr )2/3 /(r 



we can rewrite (15. 131) as 

\d 2 T + 2 e~ 6p - 2x d T - k 2 I(r) ( e - 6p - 2x sinh r) 2/3 j e = 



(5.15) 



(5.16) 



The choice of momentum scale normalization in (I5.15P is a little unusual, but convenient 
in what follows. 

In this simple one-component system, the method described in Sec. 13.11 is sufficiently 
robust For this method, first we need to find the regular behaviour at r = 0, from which 
the initial conditions will be inferred. For small r, we find to leading order q- & p~ 2x k,1/t 
and I(t) « 1, so that a regular solution of (I5.16P is found to behave as 



e reg = 1 + Y + O (r 4 ) 



(5.17) 



21 There is only one dominant and one subdominant asymptotic behaviour, and the latter is sufficiently 
suppressed for large r with respect to the former. 
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whereas the singular behaviour, which starts with 1/r, is discarded. Notice that the 
choice of momentum scale in (15 . 15[) makes the k 2 term in (15.171) independent of h(0), a 
good thing since the value of h(0) is known only from numerical evaluation of (15.111) . 

The second step of this method is to consider the asymptotic UV region. For large r 
we have, again to leading order, q~ & p- 2x 2/3 and J(r) ~ re _4r / 3 . Therefore, the generic 
dominant asymptotic behaviour is just 

c dom = 1 + (e~ 2 ^ 3 ) . (5.18) 

This tells us that we must simply search for values of k 2 for which the regular solution 
tends to zero for large r. 

We find a discrete spectrum of spin-2 states, the first of which are reported and 
compared to Krasnitz' values in Table [2] Krasnitz' results, which were obtained using a 
WKB approximation, are in good agreement with ours. As shown in Fig. (3j the spectrum 
fits nicely to a quadratic curve. A least-square fit yields 

m 2 « 0.2715 n 2 + 0.4936 n + 0.2969 . (5.19) 

As we will find in the next section, the coefficient of the n 2 term enjoys a certain degree 
of universality. 



5.3 KS 7-scalar system: Setup 

Now, we turn to the most difficult system of our paper: the seven coupled scalar fields 
in the KS background. All scalars appear to be fully coupled in the bulk, but we follow 
the insight gleaned from the "moderate UV" approximation [llj to decouple a 4 x 4 set 
of fields (the "glueball sector" ) from a 3 x 3 set (the "gluinoball sector" ) to leading order 
in the UV, as will be seen in Sec. 15.41 The system of field equations we consider follows 
from (12. 5p upon changing the radial coordinate to r. One finds 

[{d T - M) (d T -N)-k 2 e~ 2A ' 8p ] a = , (5.20) 

where we have dropped the tensor indices, and the matrices M and N are defined by 

M\ = -N a b -K\-2 e - 2x ~ 6p 51 , 

, / W a Wh\ 
r^e-Wr-^j , (5.21) 

K\ = 2e~ ip g a bc W c . 

As before, we fix the momentum scale as in ( 15. 15ft . 

Now, the matrices M and iV a priori depend on the constants P and $o- However, 
there is a linear field transformation that removes this dependence from (15.201) . This 



implies that these constants can affect the mass spectrum only by an overall change of 
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the momentum scale, which is not visible in our effective 5-d approach. Starting with 
the fluctuation vector o = (Sx, Sp, Sy, 5$, 5b, 5h\, 8h 2 ) T , the linear transformation that 
accomplishes this is a' = Ra with 

and we also rotate the matrices by, e.g., N' = RNR -1 . Henceforth, we shall consider 
the rotated matrices, dropping primes. The somewhat lengthy expressions for the rotated 
matrices K and N (that are now independent of P and $o as advertised) as well as the 
sigma-model metric G can be found in Appendix IC.ll 



5.4 KS 7-scalar system: Boundary conditions 



In this section, we will consider the asymptotic (large-r) and deep-IR (small-r) behaviour 
of the solutions of (I5.20p . which are needed to fix the boundary (initial) conditions for 
the numerical integration. 

Let us start with the large-r region. Asymptotically, the matrices and the warp term 
in (15.201) can be expanded in powers of e~ r , such thato 



K = K<® + e~ T + O (e" 2r ) , N = iV (0) + e" r iV« + O (e~ 2r ) 



(5.23) 



In what follows, we can always drop the O (e~ 2r ) terms (the reason for this will become 
clearer later on, cf. the discussion below f)5.40p ). For the sake of brevity, we write the 
matrices in block form, 

i^4x4 -?Qx3 
^3x4 -^3x3, 

and analogously for N. Then, the matrices in (I5.23P are 



K 



(5.24) 



/ 



K 



(o) 

4x4 



3(r-l/4) 


-2 1 



-"■3x3 





T-1/4 ~^ 

u 3(t-1/4) » 

1 



(5.25) 



(° 



V 2 



r-l/4 





T-1/4 / 



K 



(0) 
4x3 



K. 



(o) 

3x4 







(5.26) 
(5.27) 
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The "coefficients" may contain rational functions of r, but no exponentials. 
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and 



N, 



(0) 



4x4 



N. 



(0) 



3x3 



N, 



(0) 



4x3 




v 3x4 



A" 



(i) 

4x3 








K. 



(1) 

3x4 



-4(r - 2) 


/ 

-4(r-l) 



(i) 

4x4 



\4(r 

(i) _ 



A. 



3x3 







4(r-l) 
3(r-l/4) 




8(r-l) 
3(r-l/4) 







3(t+1/4) 
2 

9(r+l/4) 
2 

3(r+l/4) 





3(r-l/4) 





2) 



4 




1 

0/ 



It 



3(t-1/4) 




'3(r-l/4)/ 



\ 

2(r-l) 
2(r-2)y 



(5.28) 

(5.29) 
(5.30) 



(5.31) 

(5.32) 
(5.33) 



N, 



(i) 

4x3 



N. 



(i) 



3x4 



AT 



(1) 



4x4 



V 

N. 



1 


4(r-l) 


4t \ 


T+l/4 


3(t+1/4) 


3(r+l/4) 


2(t-1/4) 


4(t-1) 


4t 


3(t+1/4) 


9(t+1/4) 


9(t+1/4) 


4t 2 -5t-5/2 


16t-1 


4t 


r+1/4 


3(t+1/4) 


3(t+1/4) 








/ 


2 


8(t-1/4) 


4 


r+1/4 


t+1/4 


3(t+1/4) 


2(t-1) 


8(t-1)(t-1/4) 


4(t-1) 


t+1/4 


t+1/4 


3(t+1/4) 


2(r-2) 


8(t-2)(t-1/4) 


4(t-2) 


T + l/4 


r+1/4 


3(t+1/4) 



o \ 

-2(r-l) 
-2(r-2)/ 



(i) 

3x3 







(5.34) 



(5.35) 
(5.36) 



The transformation a — > a', cf. f !5.22p . has brought the matrices into this nice block form. 
We also need 

-2x-6p _ i j_n (5.37) 



§ + O (e--) 



as well as 



-2A- 



31/3 



i)e~-/ 3 [l + 0(e*)] 



(5.38) 



It is a useful check that the leading-order terms of these expressions coincide with the 
respective quantities evaluated in the Klebanov-Tseytlin background [24] . 

Since (15.381) is exponentially suppressed, the leading order terms of the asymptotic 
solutions will be independent of the momentum k. We note that this is just as we assumed 
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in Sec. 12.31 In contrast, in the MN system, the asymptotic behaviours of the solutions 
depend on k [cf. (14. 1 3[) . 

The asymptotic UV solutions, including the leading and some subleading terms, are 
found by iteratively solving the equations 

(d T - iV (0) ) (n) = ^ {n) + e" T iV (1) (n ~ 1} , (5.39) 

(d T - M(°') =pf T -^\ e- 2T / 3 ^ + e" T M( V n_1) , (5-40) 

where (3 = 3^ 3 k 2 /h(0), and we set <f>^ = = 0. The solutions <j>® are the leading; 

order terms of the asymptotic solutions. Again, we will drop all O (e _2T ) terms in the 
iteration. The reason for this is the following. Dropping subleading terms in the initial 
conditions leads to systematic errors in the numerical solutions, which must obviously 
be kept smaller than the relevant parts of the solutions. This implies that the initial 
conditions of all behaviours we use must be given to the same order in e~ T . We will use 
the mid-point method of Sec. 13.21 in which only the subdominant asymptotic behaviours 
are needed. Now, the leading term of the weakest subdominant behaviour goes like e~ 8r / 3 , 
whereas the strongest subdominant behaviour goes like e _T . Thus, they differ by e~ 5r ^ 3 
and we can drop O (e _2r ) corrections! 23 ! 

The somewhat lengthy expressions of the asymptotic solutions are deferred to Ap- 
pendix IC.2I 

Now, let us consider the small-r region. In order to find the independent behaviours 
for small r, we expand the matrices in fl5.20p about r = and make an ansatz of the form 

a(r) = r q (a + rai + r 2 a 2 + r 3 a 3 + ■ ■ • ) , (5.41) 

with some undetermined number q, and constant vectors a n , n = 0, 1,2, . . .. These are 
then determined recursively from the equation of motion (using computer algebra). It 
turns out that, of the 14 solutions, there exist four solutions with q = 0, one with q — 1, 
three with q = 2, four with q = —1, and one each with q = —2 and q = —3. Thus, 
there exist eight solutions, whose components are finite at r = 0, and six with singular 
component behaviour. However, the regularity condition (12.101) tells us that one of the 
solutions with q = is, in fact, singular. Therefore, we arrive at precisely seven regular and 
seven singular solutions. The independent small-r behaviours are listed in Appendix IC.3I 



5.5 KS spin-0 spectrum 

Analogously to the MN case, instead of 7 coupled 2nd order ODEs, it is more convenient 
to solve 14 coupled 1st order ODEs. Thus, let us rewrite (I5.20p as 

m) (3 • 

23 The same argument holds for the method of Sec. 13.11 where only the dominant solutions are needed. 
The weakest dominant solution goes like e~ T / 3 , the strongest one like e 4r / 3 . 
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Figure 4: det7 as a function of k 2 . For clarity, we have taken the 14th root of 
the absolute value of the determinant (leaving its sign untouched). The plot shows 
clearly that there are no zero crossings for positive k 2 . The inset shows the zeros of 
det7 for the first three spin-0 states. 



m 2 


0.185 


0.428 


0.835 


1.28 


1.63 


1.94 


2.34 




2.61 


3.32 


3.54 


4.12 


4.18 


4.43 


4.43 




5.35 


5.63 


5.63 


6.59 


6.66 


6.77 


7.14 



Table 3: Low- lying Klebanov-Strassler spin-0 mass states, extracted from zero- 
crossings in Fig. [H More values are given in Appendix |Pl 

We apply the midpoint determinant method as outlined in section 13 .21 That is, we 
compute the 14 x 14 matrix 

7 =f^ M (5.43) 

' — 'mid 

as a function of k 2 and look for zero crossings of det 7. A rough plot of det 7 is shown 
in Fig. HI We verify that there are no zero crossings for positive k 2 . The zero crossings 
of det 7 can be found by zooming in on a particular region of the plot, as shown, for 
example, in the inset. As the determinant itself changes over many orders of magnitude, 
it is useful to plot ±| det7| 1 / 14 = sign(det7)| det7| 1/,:L4 instead. The sharp turns in the 
inset are merely artifacts of this, as is apparent from Fig. El where det 7 itself is displayed 
without the 14th root, and there are no sharp turns. 

We have calculated all mass values up to m 2 = 600. The first few are listed in Table El 
for a more extended list see Table [5] in Appendix [Dl Our results exhibit two interesting 
features. First, it appears that several mass values are nearly degenerate, such as the 
two values about m 2 « 4.43 and two values about m 2 « 5.63. The origin of this in our 
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-20 



-5.64 



-5.63 



-25 



-5.7 



-5.6 



-5.5 



-5.4 



Figure 5: Here we plot det7 in the range k 2 £ [—5.8 ... — 5.3]. The inset zooms in 
on the region around k 2 = —5.63, where there are two zeros. 

calculation can be seen in Fig. [5j The det 7 curve barely seems to touch the fc 2 -axis, and 
only zooming in reveals that there are two crossings. As we do not have a dynamical 
explanation of it, it is possible that this near degeneracy is accidental, or it may indicate 
a multiplet originating from a weakly broken symmetry of the low-energy effective theory. 

The second interesting feature is the appearance, for large masses, of a periodic pattern 
of period 7 in the consecutive spectrum excitation number n. On the left hand side of 
Fig. EJ this pattern is easily visible. Hence, we split the spectrum into 7 towers of mass 
states, one of which is shown on the right hand side of Fig. [UJ From now on, we distinguish 
excitation numbers of individual towers, denoted by n t , from the excitation number n of 
the whole spin-0 spectrum. 

Like the spin-2 spectrum, the seven spin-0 spectra also exhibit quadratic dependence of 
m 2 on excitation number n t . Leaving out the first two values of each tower (i.e. nt = 1, 2) 
least square fits of the "large m 2 " behavior (using mass values up to m 2 = 600) yield 



The last digit of the leading coefficients (i.e. the coefficients of the n\ terms) can vary 
slightly if one performs the fits with more than two states from each tower left out. Still, it 
is obvious that within some uncertainties the leading coefficients all coincide rather well, 




0.271 n 2 + 0.774 n t + 0.562 , 
0.270n 2 + 0.928 n 4 - 0.430 , 
0.275n? + 0.769/it + 1.921 , 
0.272n 2 + 1.017^ + 1.023 , 
0.272 n 2 + 1.119 n t + 0.337 , 
0.271 n 2 + 1.150/14 + 0.648 , 
0.273 n 2 + 1.082 n t + 2.544 . 



(5.44) 
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Figure 6: On the left, we see that for high m 2 , it is easy to identify the periodicity 
"empirically" as 7. On the right, we have extracted one tower (i.e., every 7th point) 
from the full spectrum, enumerated states in that tower by a parameter nt (the 
"excitation number") and fit those points to a quadratic curve starting with the 
third point. As shown in the inset, m 2 ~ n 2 is not a good description for very low 
n t , unlike in the spin-2 case. 

both amongst each other and also with the coefficient of the spin-2 tower (15.191) . leading 
to roughly universal behaviour at large nt- 

5.6 KS 7-scalar system: Quadratic confinement 

If one plots the experimental values for squared meson masses m? (e.g. for the first few 
resonances of the p meson) against their consecutive number n, the data obeys m 2 ~ n 
to good accuracy] 24 ! One might call this "linear confinement". Karch, Katz, Son and 
Stephanov (KKSS) [T7] emphasized that there are no known supergravity backgrounds 
that reproduce this experimental fact about mesons from AdS/CFT. 

In fact, the claim of KKSS was even stronger (p. 13): that for all known supergravity 
backgrounds m 2 ~ n 2 , which is "quadratic confinement". Extending these arguments to 
our setting, one might expect to see quadratic confinement in, say, the Klebanov-Strassler 
background. A priori, given the complexity of fluctuations around the Klebanov-Strassler 
solution, one could have expected that no such general statement could be made except 
possibly in some extreme asymptotic regime. 

This expectation turns out to be wrong. Although we cannot be absolutely sure from 
our analysis, the KKSS claim that m 2 ~ n 2 seems resoundingly apparent in Fig. [U] and 
the fits (I5.44p . Since the KS theory is in any case not QCD, linear confinement is not 
crucial to have, but our results strengthen the case for finding controllable models that 
do exhibit linear confinement (the metric-scalar model given in [T7j seems promising, but 

24 Since there should be no confusion here, we omit the subscript "t" on n t from now on. 
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so far has not been shown to be the solution of any known theory). 

Of course, one would expect that much of the low-energy dynamics of the theory would 
be determined by the lowest-lying mass states, and they have a richer structure than the 
large-n asymptotics — a glimpse of this can be seen in the inset on the right hand side of 
Fig. El Since the KS theory incorporates spontaneous breaking of chiral symmetry, this 
should have left its mark on our spectrum. We do observe some interesting patterns, but 
we have not been able to link them conclusively to symmetry breaking. This deserves 
further study. 

As alluded to in the introduction, operator mixing is even more subtle here than in 
standard AdS/CFT. We have computed the spin-0 mass spectrum, but we have not said 
what the mass states are composed of in terms of dual gauge theory operators. Some of 
the problems of identifying the mass states were already addressed in section 12.31 i.e. the 
ambiguity in the normalization of the dominant asymtptotic solutions and the freedom 
to add to a given dominant solution multiples of asymptotic solutions of equal or weaker 
strength. Thus we do not attempt to give a further interpretation of the mass states. We 
just content ourselves by noticing that there is actually a quantity that is easy to obtain 
with our method, namely the coefficients from (j2.28j) . These can be determined 
straightforwardly by solving f)3.8p at the various mass values. This gives information on 
the composition of the "cavity mode" (that we called <X\ in chapter [2]) corresponding to a 
mass state in terms of our basis of subdominant solutions. However, as the interpretation 
in terms of operator mixing is not clear to us, we refrain from determining these coefficients 
explicitly here. Moreover, also the quantities Z\ ti , defined in (12.331) . could be determined 
in principle. They correspond to the residues of the mass poles if N\, introduced in (I2.32|) . 
is really one, as in the case of asymptotic AdS spaces (see Appendix |A"|) . 

5.7 Comparison to hardwall approximations 

In this section, we would like to give some arguments supporting our conviction stated in 
the introduction that the detailed form of the spectrum depends crucially on the details 
of the IR dynamics in the gauge theory. Holographically, this means that it is important 
to use the full Klebanov-Strassler background to determine the low lying mass eigenval- 
ues. To illustrate this, we compare our spectrum with a naive toy hardwall model. In 
particular, one might hope to find the right spectrum by using the Klebanov-Tseytlin 
approximation of the KS background. This is roughly equivalent to cutting off the KS 
background at some value Tir of order one (the two solutions differ strongly only for 
r < 0.8). In order to use the determinant method described in Sec. 13. 1[ one would like to 
impose seven independent initial conditions at 7ir and then solve the equations numeri- 
cally from the IR towards the UV. One again demands that there is a linear combination 
of those numerical solutions that does not contain any dominant solutions adom when 
expanded in a basis of asymptotic solutions. Lacking any guiding principle for choos- 
ing these initial conditions, we arbitrarily impose that all components of the seven scalar 
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values from Table [3] 


T ra = 0.5 


T m = 0.75 


Ttr = 1 


0.19 


0.27 


0.30 


0.34 


0.42 


0.56 


0.63 


0.71 


0.83 


0.86 


1.03 


1.14 


1.29 


1.39 


1.62 


1.83 


1.63 


1.78 


1.89 


2.03 


1.94 


2.32 


2.42 


2.57 


2.35 


2.43 


2.56 


2.75 


2.61 


3.26 


3.39 


3.70 


3.33 


3.58 


3.80 


4.09 


3.53 


4.06 


4.22 


4.40 


4.13 


4.30 


4.60 




4.17 


4.54 


4.81 




4.43 (2x) 


4.71 







Table 4: Comparison of our mass values with a toy hardwall model as described 
in the text (only values up to m? = 5 are included). This demonstrates that the 
spectrum indeed strongly depends on the IR dynamics. 

fields should vanish at tir, whereas for each initial condition one of the first derivatives are 
taken to be unity. This gives seven independent initial conditions. The resulting spectrum 
is contained in Table H] for several values of tir, together with the correct mass values. 
Obviously, the "toy spectra" depend on the value of tir, but they also differ strongly from 
the correct values. We believe that this holds more generally, even if one had a different 
set of initial conditions than these, whether or not they are physically motivated. 

To summarize, we do not see any way of implementing a hard-wall approximation that 
reproduces the correct spectrum of Table [HJ 

6 Outlook 

We have computed spin-0 and spin-2 mass spectra for the Maldacena-Nunez (wrapped D5- 
brane) and Klebanov-Strassler (warped deformed conifold) backgrounds. Although fairly 
complicated, there are some simple features. For example, for large excitation number n, 
the spin-0 states in the KS theory organize into 7 towers, each of which displays quadratic 
confinement, in agreement with the claims of [17]. ("Large" is not terribly large; around 
n > 3.) For low excitation number n, there is a rich structure of mass values, which 
exhibits interesting patterns of near degeneracy. 

For the Maldacena-Nunez background, our numerical results confirm the rough fea- 
tures of the analytic spectra we derived in [TT] and improved upon here; the discrete mass 
spectrum has an upper bound, which is quite different from the KS case. 
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The situation with respect to other studies of glueball spectra in these backgrounds, 
based on various truncations, has not changed from [IT] . The one result we have been 
able to compare our results to is the set of mass states for spin-2 glueballs in KS, where 
the first few values were computed in the WKB approximation by Krasnitz in [40j, and 
are roughly correct, as in Table [2] 

It would be very important to understand the details of mixing in this context; without 
further information along the lines we have suggested, we cannot determine the composi- 
tion of the mass states we find. On the good side, the spectrum we have found is unique, 
in that the variables used to find it are gauge- invariant, and we impose the boundary 
conditions one really wants to impose, rather than approximations thereof. 

These results and methods could straightforwardly be applied to other gauge/gravity 
dual pairs. Apart from other J\f = 1 dualities, it would be interesting to apply these 
methods also in Af = gauge theory. There, the existence of a simple superpotential W, 
which greatly simplified our task, is not guaranteed but occurs e.g. for Af = backgrounds 
with "pseudo-Killing" spinors (cf. for example [21] ). Here is a concrete question: do the 
meson spectra computed in other backgrounds contain glueball admixtures from the 
in the DBI action? 

Independently of supersymmetry breaking, including probe branes in Klebanov-Strassler 
a la [18] [12] would be interesting to study using our methods. 

Finally, the greatest challenge for holography in presently known confining back- 
grounds is to compute correlators from the Klebanov-Strassler background. Some progress 
was made in [3 HTJ HQ], H2} H3]. We considered the Krasnitz (extreme high-energy) limit 
in [11] and found some analytical solutions. However, those solutions were checked to be 
valid around k 2 ~ 10 6 , whereas in this paper, we were interested in IR physics, and thus 
the Krasnitz approximation was not available to us. In fact, we only considered solutions 
up to k 2 ~ 2000, and going much higher presents a numerical difficulty, as one needs to go 
increasingly far in the UV for the solutions to become asymptotic, so that the determinant 
method is applicable. This is a superable challenge for connecting the high-energy and 
low-energy regimes, which might ultimately be the most useful aspect of holography for 
confining gauge theories. We hope to return to these issues in future publications. 

Acknowledgements 

It is a pleasure to thank Massimo Bianchi, Andreas Karch, Emanuel Katz, Albion Lawrence, 
Scott Noble, Carlos Nunez, Henning Samtleben, David Tong and Amos Yarom for helpful 
discussions and comments. This work is supported in part by the European Community's 
Human Potential Program under contract MRTN-CT-2004-005104 'Constituents, funda- 
mental forces and symmetries of the universe'. This research was supported in part by 
the National Science Foundation under Grant No. PHY99-07949. The work of M. B. is 



32 



supported by European Community's Human Potential Program under contract MRTN- 
CT-2004-512194, 'The European Superstring Theory Network'. He would like to thank 
the Galileo Galilei Institute in Florence for hospitality. The work of M. H. is supported 
by the German Research Foundation (DFG) within the Emmy Noether-Program (grant 
number: HA 3448/3-1). Both M. B. and M. H would like to thank the KITP in Santa 
Barbara for hospitality during the program "String Phenomenology" . The work of W. M. 
is supported in part by the Italian Ministry of Education and Research (MIUR), project 
2005-023102. 

A Some 2-point functions in AdS/CFT 

In this section, we review how the 2-point functions for some known cases of AdS/CFT 
arise as sums over the spectrum of bulk eigenf unctions, using expressions developed in 
the main text. 

A.l AdS bulk 

Let us start with a set of n free massive scalar fields on (d + l)-dimensional AdS bulk 
space. Using the language of Sec. [2} we will derive the 2-point function as a "sum" over 
the spectrum of bulk eigenfunctions. (For pure AdS the spectrum is continuous, so the 
"sum" is in fact an integral. This will be different in the next subsection.) When the 
smoke has cleared (in ( 1A.14|) ). the connection to the standard computation in terms of 
the IR-regular Bessel function K u (kz) will be evident. 

The bulk equation of motion (12. 5ft (in d-momentum space) is 

(di-^d z -^- k *y=o, (A.i) 

where z = e~ r is the radial coordinate. We have omitted the field indices, and m 2 is to be 
understood as a diagonal matrix with entries mf , i — 1, . . . , n. The conformal dimensions 
of the operators Oi dual to the components of <fi are related to the mass parameters rrii 

bjH 

d [d? 
Ai=- + a f , ai = \j- + m 2 . (A.2) 

Small values of z give the asymptotic UV region, whereas large z describe the bulk interior 
(IR region). 

25 In general one needs to consider = | ± on, but for the purposes of this discussion we restrict to 
the upper sign. For simplicity, we also consider generic (non-integer) values of ct^. 
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Conventionally normalized asymptotic solutions of ( 1 A. II) ar 



(f>^(z,k) = SfT{l-a i )(^\ * z d / 2 \_ a Xkz) = 5*z d/2 - ai + -- - , (A.3) 

<M(z,k) = str(i + oi)(^ \ d / 2 \ a Xkz) = 5U dl2+a * + ■■■ , (A.4) 

where the l a are modified Bessel functions, and after the second equal signs we indicated 
just the respective leading behaviours. The powers of k in front of the solutions are nec- 
essary in order to make the leading terms k- independent, while the remaining coefficients 
are conventional normalizations. 

For asymptotically AdS bulk backgrounds, (lA.ip still holds in the asymptotic region. 
Therefore, the leading behaviours of the asymptotic solutions remain as in (jA.30 and 

dA3D, 

0dam?(*> k) = 51 z d ' 2 -^ + ■■■ , cj> snh «(z, k) = 51 z d ' 2+ai + ■■■ . (A.5) 

The general form of the matrix 112. 30[) in asymptotically AdS spaces follows easily. 
One finds 

Z i:j = 2oii 5 i:j . (A. 6) 

A general result of holographic renormalization [MJ ESI EE] is that (the non-local part 
of) the exact one-point function (Ci) exact is 

{O i ) tm * = 2a i d i ^ (2A, - d) * . (A.7) 
Comparing this with 1 12.241) . one finds 

Y ij = 2oL i S ij = Z i j , (A.8) 

which, together with (12.32ft . implies N\ = 1. 

Returning to the pure AdS bulk, ( 1A.1I) admits a continuous spectrum of regular and 
subdominant solutions for k 2 = —A 2 , A > 0. The eigenfunctions are given bvF^l 

^) = ^V\z d / 2 J a X\z) (A.9) 

and satisfy the orthogonality relation 

CO 

—z-WfaW ■ <p xtj (z) = 5(X - \')5 i3 . (A.10) 
o 

Considering the small- z behaviour of ( 1A.9I) . one can read off the response coefficients d\i, 



' \ \ on /T 

J \2 J 1(1 • n,) V 



26 We have reinstated the field index a. 

27 The generic label A for the eigenfunctions used in the main text [c.f. (|2.14[) ] is replaced here by two 
indices, A = — k 2 and i = 1, . . . ,n. As before, the upper index a is the vector component index. 
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Hence, after using (12.331) and N\ = 1, one obtains for (12.341) 



oc 



2 2(i-a,) f A 2<,,+i r(l-a ( ) /kV' 



(OmOA-k)} = 4 jTRr y dA— = -« B 2-^ (A.12) 



Notice that the second equality holds only after analytic continuation, because the integral 
does not exist if aj > 0. This is equivalent to adding an infinite contact term to the integral 
over the spectrum. For example, if < «j < 1, we rewrite the integrand as 

X^-k 2 ^— (A.13) 



k 2 + A 2 fc 2 + A 2 

and add a counter term that cancels the first term on the right hand side. 

The finite result flAT2l can also be obtained in the usual way, i.e., by considering a 
regular solution of (1A.1I) . 



CO 



2d f k\ ai ^O l ~ a i r \ ai+1 

= *F(5J (2) = s " d fW)J AX ^ z " 2 ^ {Xz) ' {AM) 



reading off the relation between the response and source coefficients, d and d, with the 
help of 

KJx) = 2 a - 1 T(a)x- a (l + ...)- 2~ a - l T ( 1 ~ a K a (l + . . .) (A.15) 

a 

and then using ( 1A.7I) . It is interesting to note that the small-z behaviour of the inte- 
grand in (IA.14I) is only subdominant, while the expression in terms of K a . contains also 
a dominant piece. The explanation of this apparent discrepancy is the UV part of the 
spectrum (large A), i.e., one cannot find some sufficiently small z such that Xz is also small 
for all values of A. Hence, loosely, the UV-part of the spectrum generates the dominant 
behaviour. 



A. 2 Active scalar in GPPZ flow 

As a second example, let us take the active scalar in the GPPZ flow [H]. Its gauge 
invariant bulk equation of motion reads 



u(l - u)d 2 u + (3« - 2)d u + - 1 + T 



, (A.16) 



where the radial coordinate u is defined by u = 1 — e 2r , and the warp factor is e 2A = 
u/{l-u). 

Eq. (jA.160 admits a discrete spectrum of regular and subdominant eigenfunctions, 
with mass squares 

m 2 n = 4n(n + l), n = 1,2,3,.... (A.17) 
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The normalized eigenfunctions are 



'2(2n + 1) 
n(n + 1) 



(1 



u 



,3/2. 



du 



P n {2u - 1) 



(A.18) 



where P n are Legendre polynomials! I One easily finds the response coefficients 



d n = y/2n(n + l)(2n + 1) . 



(A.19) 



Thus, using a = 1 (we have a A = 3 operator) and (1A.8D (which holds for asymptotically 
AdS spaces), we obtain for the 2-point function (12.341) 



(0(k)0(—k)) = J2 



ln(n + l)(2n + 1) 
^ k 2 + 4n(n + l) 

n=l v ' 



+ C.t. . 



(A.20) 



Clearly, the sum in (lA~20j) does not converge, so that there are again infinite contact terms. 
It is instructive to compare ()A.20j) with the finite result from holographic renormalization 

sua, 



(0(k)0(-k)) 



■0 



3 + VT^¥ 



if, 



V>(l)-V(2) 



(A.21) 



where ip{z) = [\nT(z)}'. Using the identity 

°° / 1 

if>(x) -i/>(y) = ( — 7 

n=o \y 



we obtain from (1A.211) 

(o(*)f(-*)) = -E 



2n + 1 



+ n x + n 



(A.22) 



2 ^n(7i + l)[A; 2 + 4n(n + l)] 
n(n+l)(2n+l) 



n=l 



k 2 + 4n(n + 1) 



2(2n+ 1) + t 



2n + l 
2n(n + 1) 



(A.23) 



Again, we see that holographic renormalization makes a precise choice for the contact 
terms and, therefore, picks a renormalization scheme and scale. Indeed, the finite result 
( 1A.21I) satisfies the renormalization conditions 



(0(k)0(-k))\ k2 __ 







dk 2 



(0(k)0(-k))\ k2 __ 



. 



(A.24) 



2<x 



As in [221 HQ], regular and subdominant solutions of (|A. 16|) are given by Jacobi polynomials P n l 1 ( 
which are proportional to dP„(z)/dz. 
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B Matrices for the MN background 



The entries of the matrices M and N of Eq. (14.71) for the regular MN background (c = 0) 
are 

Mn = -[(8p 2 - 4p+ 1) + e" 4p (-64p 3 + 112p 2 - 48p + 2) + 
+ e" 8p (256p 4 - 896p 3 + 248p 2 + 108p - 17) + 
+ e" 12p (1536p 4 - 736p 2 + 28) + 
+ e~ 16p (256p 4 + 896p 3 + 248p 2 - 108p - 17) + 
+ e" 20p (64p 3 + 112p 2 + 48p + 2) + e" 24p (8p 2 + 4p + 1)]/ 
/{p(l - e- 4 ') 4 [l - e" 4 "(4p + l)](4p - 1 + e" 4 ")} 
Mia = [(3p - 1) + e" 4p (16p 3 - 48p 2 + 12 p + 2) + e~ 8p (96p 3 - 30p) + 
+ e" 12p (16p 3 + 48p 2 + 12p - 2) + e~ 1(5p {3p + 1)] x 
x 4e- 2 " /{p(l - e- 4p ) 3 [l - e~ 4p (4p + l)](4p - 1 + e" 4 ")} 
M 2i = -[1 - e- 4 "(4p + l)](4p - 1 + e- 4p )[{p - 1) + e" 4 "6p + e- 8 "(p + l)]x (B.l) 

x 4e^/[p(l -e- 4p ) 5 ] 
M 22 = -[(2p - 1) + e" 4p (32p 3 - 128p 2 + 56p - 2) + 

+ e" 8p (-256p 4 + 960p 3 - 256p 2 - 118p + 17) + 
+ e" 12p (-1536p 4 + 768p 2 - 28) + 
+ e" 16p (-256p 4 - 960p 3 - 256p 2 + 118p + 17) + 
+ e- 20 "(-32p 3 - 128p 2 - 56p - 2) + e- 24p (-2p - 1)]/ 
/{p(l - e-^) 4 [l - e^(4p + l)](4p - 1 + e"*)} 
M 33 = -8p(l - e" 4 2 /{[l - e" 4 "(4p + l)](4p - 1 + e" 4 ")} 
Mia = M 23 = M 3 i = M 32 = , 

and 

Nu = -[1 + e~ 4p (16p 2 - 32p + 4) + e~ 8p (96p 2 - 10) + e" 12p (16p 2 + 32p + 4) + 

+ e- 16 "]/[p(l-e- 4 Y] 
N 12 = [(p - 1) + e~ 4 " 6p + e- 8 "(p + 1)] x 4 e~ 2f> /[p(l - e" 4 ") 3 ] 
A^ 2 i = -[(3p - 1) + e- 4p (16p 3 - 48p 2 + 12p + 2) + e- 8p (96p 3 - 30p) + 

+ e~ 12 "(16p 3 + 48p 2 + 12p - 2) + e- 16 ^(3p + 1)] x 4 e" 2 ^ /[p(l - e" 4 ") 5 ] (B ' 2) 
N 22 = -[(2p - 1) + e" 4p (-16p 2 + 28p - 4) + e" 8p (-96p 2 + 10) + 

+ e- 12 "(-16p 2 - 28p - 4) + e- wp (-2p - l)]/[p(l - e" 4 ") 4 ] 
iVis = N 23 = N 31 = N 32 = N 33 = . 
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C Bulky KS stuff 



C.l KS matrices 

We present here the explicit expressions for the 7x7 matrices appearing in Section 15.31 
To shorten the formulae, a number of abbreviations will be used. First, 

c = coshy = cothr , s = sinh?/ = — (sinhr) -1 , (C.l) 

where y denotes the background field of Section 15.11 Second, we introduce 

B 1 = tc-1, B 2 = rs 2 - c , (C.2) 



and 



A x = h(r) (AsB 2 y 1/3 = h(r) sinh r (2 sinh 2r - 4r)~ 1/3 , (C.3) 
A 2 = -Ax (cB 2 - - l -sB x B 2 . (C.4) 

Let us consider the behaviour for small and large r of Ax and A 2 . As h(0) is a finite, 
positive constant, one obtains 

Ax(0) = ±pWh{0) , A 2 (0) = ^(0) . (C.5) 



For large r, starting from 



(C.6) 



one obtains 

Mr) « ~ e- (r - ~) , A 2 (r) * ~ [r + i) . (C.7) 

With the abbreviations (ICHf - rfCliL the (rotated) matrices K\ = 2e~ ip g a bc W c and 
N a b are given by 

/ n — - — o o s2Bi s2t 

1 u u 2A X B 2 U U 2A X B 2 2A 1 B 2 



2(l + 2c£ 2 ) 2 (^+cAiB 2 ) 2c5 2 + l 2s(t + 2B 2 ) -2s 



K 



AiB 2 

-rV o 



„2, 



s 



ooo Ul , i 

2sB 1 " -sBx 2(A2 + C ^ li?2) ' 

\-2s{t + 2B 2 ) -2s -s(r + 2£? 2 ) -2(1 + 2cB 2 ) V ^2+£pM J 

(C.8) 
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/ 



N 



2c(A 2 +cA 1 B 2 ) 
A 2 
_2cA 1 _ 
3A 2 

-2(l+2cB 2 )(A2+cA 1 B 2 ) 
A 2 



2s(A 2 +cA 1 B 2 ) 
A 2 

2sB 1 (A 2 +cA 1 B 2 ) 
A 2 

2s{t+2B 2 )(A 2 +cA 1 B 2 ) 
A 2 





cs 




2A 2 






A 2 


&A 2 B 2 


4Ai(l+2cB 2 ) 


s(l+2cB 2 ) 


A 2 


2A 2 








isA\ 


s 2 




9 4„ 


4sAiBi 




A 2 


2A2 


AsA 1 {t+2B 2 ) 


s 2 (r+2B 2 ) 


A 2 


2A 2 






-(l + 2cB 2 ) 



s(r + 25 2 ) 



s(A 2 +cAiB 2 ) 
A 2 

3A 2 

2sA 2 (t+2B 2 )+sA 1 B 2 (1+2cB 2 ) 
A 2 



s 2 A 1 B 2 _ 
A 2 C 
s 2 A 1 B 1 B 2 

A 2 

2A 2 (1+2cB 2 )+s 2 A 1 B2(t+2B 2 ) 
A 2 



_csf_Bi 

2A 2 
s 2 B 1 
6A 2 B 2 
s 2 B 1 (l+2cB 2 )+4csA 2 
2A 2 





CS 2 T \ 
2A 2 \ 



2A 2 

s :i B 2 



2A 2 

s 3 B 1 (t+2B 2 )+2A 2 (2s 2 +1) 
2A 2 



6A 2 B 2 
s 2 t(1+2cB 2 ) 
2A 2 



2A 2 
s i rB 1 1 1 
2A 2 ^ L 

s 3 t(t+2B 2 ) 

2A 2 



(C.9) 



The matrix M is given by 



M = —N — K + 
where I denotes the 7x7 unit matrix. 



(CIO) 



Finally, we also need the sigma-model metric for the rotated fluctuation fields. It 
transforms as G' = (i? _1 ) T Gi? _1 , where R is the linear transformation matrix that leads 
to ( 15.221) . and the superscript T denotes the transpose. Explicitly, we find 



G' 



(\ 

6 




\0 













iP 2 e* " 2:r cosh(2y) 








lp2 e <fo-2x sinh ^ 22/ ) 



1 

\ 
















1 p2 e $ -2x 







sinh(2y) 



(C.ll) 



i P 2 e <Ko-2x coshpy) J 



where for x and y one should substitute the respective backgound solutions. 
C.2 KS 7-scalar system: Asymptotic solutions 

In this appendix, we list as reference the asymptotic solutions of KS the 7-scalar system. 
They are found by iteratively solving (I5.39P and (15.401) . 



For convenience, we abbreviate 







31/3 

Mo)' 



(C.12) 



Moreover, O3 and O4 denote 3 and 4 zero vector entries, respectively. 
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The dominant asymptotic (large r) solutions, up to and including terms of order e T//3 , 



are 



At + 1 



/-12\ 
4 
12 


V°3/ 



9(3 



32(4r + 1) 



,2t/3 



/ 6(5 + 4r) \ 
-(9 + 4r) 
-6(5 + 4r) 


\ o 3 J 



24 



4r+ 1 



e -/3 



/ 4 \ 

1 

T- 1 

v 2 -v 



+ 



27/3 2 



256(4r + 1) 



/-24r 2 - 48r - 63/2\ 
8r + 9 
24r 2 + 48r + 63/2 


o 3 



/o 4 \ 



1 

V 1 / 



+ 32 



V 

/ o 4 \ 

2 

2 - 2r 
V-l - 2rJ 



27P(Ar + 5) /3 
8(4r+ 1) 



/ 4 \ 
1 

r- 1 
V2-r/ 



(C.13) 



/ 2 \ 

-2/3 
4r - 2 


v °3 y 



256 



/ 4 \ 

8r 2 - 30r + 45 
-6r 2 - 39r + 243/2 
y 6r 2 + 6r - 81 J 



4r + 1 



/4r+ 13 \ 
2r - 7/2 
12r-9 


V °3 J 



L 

32 



/ 36r + 63 \ 
St- 18 
24r 2 + 12r- 279/2 
72r + 42 
3 



-T/3 



4r + 1 



4r + 1 



-t/3 



/ 

/ 1/2 \ 

-1/6 

T - 1 

1 

V °3 / 

/ V2 \ 

-1/6 
r- 1/4 


v °3 y 

/ o 4 \ 

2r+ 1 
3r + 3/2 

V 9/4 y 



o 4 

64r 2 - 104r - 6 
120r 2 - 246r - 99 
^-24r 2 + 174r + 99y 



(C.14) 
\ 



(C.15) 
(C.16) 



(C.17) 



(C.18) 
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and 



-r/3 



/0 4 \ 

4 
9 

v-v 



(C.19) 



The subdominant asymptotic (large r) solutions, up to and including terms of order 



-81-/3 



(why exactly e 8r ' 3 is explained in section l5~4"l) . are 





/ 


4 \ 




/ 


4 \ 


e" T 




1 






4r + 11 






T 


16 




2r 2 + r - 12 




V 1 


-v 




V 


-2r 2 + r + 85/4y 



-2- 



+ 



20 



/ 32r-29 \ 
|(16r + 37) 
-7(16r-7) 


3 



+ 



27{3 2 
6400 



-7r/3 



0, 



\ 



+ 







16000(4r+ 1) 



-8r/3 



300r 2 + 1625r + 627 
80r 3 - 33r 2 - 1131r - 13649/10 
^-^(1600x3 + 2440r 2 - 35920r - 27383) J 

I 3(35200r 3 - 27760r 2 - 316608r + 42739) \ 
4(800r 3 + 4240r 2 - 15228r - 33961) 
-3(169600r 3 + 219680r 2 - 705996r - 67973) 



3 



(C.20) 



/0 4 \ 


1 

v-v 



/ 



9(3 



-5r/3 



\ 



4 
1 

r+ 1/ 
\-r + 5/8j 

( 



+ e 



-It 



+ 



3^ 
640 



-7r/3 



/ -3 \ 

-13/6 
1 


V °3 J 
4 

195r + 2989/12 
^(3150r 2 + 1505r- 5591) 



\ 



OU v 

\ v - T ^(11700r 2 - 5260r 



+ 



3/3 

400(4r + 1) 



-8r/3 



28993) J 

( -3(40r 2 + 682r + 69) \ 
-(4r+ ll/2)(140r + 57) 
-3(620r 2 - 709r- 117) 


3 



V 



(C.21) 
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V 



1 

-1 

6r-3 
-4r + 9 
3 



7 



+ 



(3 2 



3(3 _ 
800 



2r 



V 



204800 



-8t/3 



/ -80r 2 - 32r - 291 \ 
-40r 2 - 356r - 5807/6 
1680r 2 + 2732r - 1634 
-50(4r- 15)(4r + 9) 

3 J 

( o 4 \ 

2000r/3 
400r 2 + 260r - 547 
\-A00t 2 + 2260r - 297 J 

( -34560r 3 - 111024r 2 + 259556r - 497927/2 \ 
11520r 3 - 104688r 2 - 1489028r/3 - 2917225/6 
470880r 3 + 1839672r 2 - 135995r - 766688 
36000r 2 - 1110600r - 1043879) 
3 



+ _^_ e -7r/3 

250 



-3(57600r 3 



/ 



e~ 4r / 3 
At + 1 



3 

-1 
12r 
-4(4r + 1) 
3 



+ 



3(3 



,-2r 



20(4r+ 1) 



25(4r+ 1) 



+ 



9f3 2 



10240(4r + 1) 



-8t/3 



/ -(12r 2 -51r-23) \ 
-i(144r 2 -56r + 53) 
2(56r 2 + 17r-4) 
-5(4r+ l) 2 
V 3 J 

( 4 \ 

-50 

-80r 2 - 74r + 49 
^-SOr 2 + 226r - 51 J 

( -1152r 3 + 2720r 2 + 7122r + 1763 \ 
-384r 3 - 992r 2 /3 + 650r + 361/3 
256r 3 + 10444r 2 - 1392r - 4169/4 
-3(320r 2 + 200r - 343) (4r + 1) 
3 



-7r/3 



-2t 



At + 1 



/ At + 1/5 \ 
2r + 23/30 
-At - 1/5 


3 



+ 



3(3 



160(4r+ 1) 
/ 

„-7r/3 



-8t/3 



/ 80r 2 + 144r + 5 \ 
|(20r 2 + 36r + 11) 



-(80r 2 + 144r + 5) 


V °3 J 

4 \ 

1/2 
•!(5r + 4) 



V-tm( 10 ^- 17 )/ 
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and 



-8t/3 



30(4r + 1) 



/ 3(160r 2 - 172r + 1) \ 
-(160r 2 + 308r+ 121) 
-6(260r 2 - 107r - 16) 
-450(4r + 1) 
3 



(C.26) 



C.3 KS 7-scalar system: Small-r behaviour 



In this appendix, we provide the independent small-r behaviours. In the calculations, 
only the regular solutions are needed. We include subleading terms up to and including 
order t 2 , as they are needed to avoid systematic errors in the initial conditions. For 
completeness, we also list the singular solutions. In the following, we abbreviate 



A 1 = A 1 (0) = -3 1 ' 3 h(0) . 



(C.27) 



The seven regular small-r solutions are 
/ 4 \ / 2 \ 



4 
1 


-4 
-2 
-1QA 1 




+ O (r 4 ) , r 2 



2 
1 

(4/3)r 


-2 


\(4/3)r/ 



+ 0(r 4 ) , r 2 



13 

3/2 
(32Ai/3)r 
20 
17 




+ 0(r 4 ) , 



M 




/ -3A; 2 /20 + 9/(80Ai) \ 




M 




/-A; 2 /60 + l/(120Ai) \ 









-3^/40 + 9/(160^) 









-Jfe 2 /120 + l/(240Ai) 




l 






















+r 2 







1 


+r 2 


A; 2 /6 + l/(6Ax) 









3A; 2 /20 + ll/(80Ai) 









fc 2 /60 + 19/(120Ai) 
























w 




V (V6)r ) 




w 




V o j 





i 
o 
o 




V o / 



( -2P/5 - 3/(8A) \ 
13fc 2 /60 + l/(16Ai) 



-P/10 + 3/(8A 1 ) 


V o / 



+ 0(r 4 ) , 
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/1\ 




2 


W 



( k 2 /5- l/10 + 31/(120Ai) \ 
k 2 /60 - 1/20 + 31/(24(L4i) 


2/(3^) 
3P/10 - 9/10 + 109/(120Ai) 





+ O (r 4 ) . 



(C.28) 



The seven singular small-r solutions are 



/ 



\ 






8A1/3 





1 



+ 0(t) 



1 

r 



J 



n\ 




2 
2 


\0J 



+ 0(t) 
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T 
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12 
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0(t) 
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T 
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( 1 \ 
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-1 


V / 






1 






0{r) , -3 





1 




-T 

V 1 / 



+ 



/ 


-4^73 + 4^/9 + 1/6 








/ \ 



12 
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AA X 
\ / 
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O(r) 



+ 0(t) 



-AAik 2 









22Ai/5 




2/3 
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/ 



16Ai/c 2 /3 + 20Ai/3 




V 



\ 



2/9 



0(r) 



(C.29) 



D KS spin-0 spectrum 

Table [5] gives a somewhat more extensive collection of the spectrum we have found (we 
computed the spectrum up to m 2 = 600). 
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n 


2 

777 


n 


2 

777 


n 


2 

777 


77, 


2 

777 


n 


2 

777 


1 


0.185 


16 


5.63 


31 


12.09 


46 


21.33 


61 


32.30 


2 


0.428 


17 


5.63 


32 


12.99 


47 


21.58 


62 


33.04 


3 


0.835 


18 


6.59 


33 


13.02 


48 


22.10 


63 


34.82 


4 


1.28 


19 


6.66 


34 


13.31 


49 


23.53 


64 


35.21 


5 


1.63 


20 


6.77 


35 


14.23 


50 


23.95 


65 


35.54 


6 


1.94 


21 


7.14 


36 


15.03 


51 


24.24 


66 


37.65 


7 


2.34 


22 


8.08 


37 


15.09 


52 


25.94 


67 


38.17 


8 


2.61 


23 


8.25 


38 


16.16 


53 


26.32 


68 


38.47 


9 


3.32 


24 


8.57 


39 


16.89 


54 


26.67 


69 


39.32 


10 


3.54 


25 


9.54 


40 


17.03 


55 


27.30 


70 


41.15 


11 


4.12 


26 


9.62 


41 


17.44 


56 


28.95 


71 


41.79 


12 


4.18 


27 


9.72 


42 


18.61 


57 


29.25 


72 


42.01 


13 


4.43 


28 


10.40 


43 


19.22 


58 


29.62 


73 


44.22 


14 


4.43 


29 


11.32 


44 


19.40 


59 


31.57 


74 


45.01 


15 


5.36 


30 


11.38 


45 


20.79 


60 


31.93 


75 


45.19 



Table 5: Klebanov-Strassler spin-0 spectrum, first 75 values. 

E Numerics 

In this section, we give a few details on how we computed the KS mass spectrum and 
how sensitive this spectrum is to numerical uncertainties. (Not very sensitive at all, as 
we now argue.) A quick summary is given in Table El The main program parameters are 
r iR; r uv, r mid and n imposed . Let us begin by explaining what they are. 

The first two are self-explanatory: the actual radial interval of the Klebanov-Strassler 
background is r = [0,oo], but to put the equations on a computer we need to impose 
some cutoffs@ Thus, tir, is an IR cutoff and tuv is a UV cutoff. The parameter r mi d is 
the midpoint for the midpoint determinant method (see section I3T2I) . i.e. the position in 
r where we compute the matrix 7^. We will explain the parameter ni mpose d in a moment. 

To integrate the 14 1st order ODEs between tir and Tyy, starting from analytically 
known boundary conditions given in previous appendices, we use a standard adaptive 
Runge-Kutta 4th-5th order solved, implemented in Maple, Fortran and C. For purposes 
of checking indiviual mass values, the Maple code is sufficient. For larger scans of param- 
eter space, the Fortran and C codes were useful. For a typical run of 1, 000 values of k 2 

29 0ne could map the domain to e.g. [0,1], but in simple examples we did not see any improvement 
by doing so. Also, it might seem confusing that we need to impose an IR cutoff when the background 
is smooth in the IR. Just as in AdS, though, the equations of motion are singular for r = 0, but the 
singularity is regular. However, some of the regular boundary conditions coincide for r exactly zero, so 
to implement them faithfully we need tir > 0. 
see e.g. www.nr.com 



45 



TlR 






^imposed 


2 

m 


0.001 


20 


1 


500 


0.18500 


0.2 


20 


1 


500 


0.18501 


0.001 


11 


1 


500 


0.18500 


0.001 


25 


5 


500 


0.18500 


0.001 


20 


1 


1 


0.1850x 


0.001 


20 


1 


500 


1.28107 


0.2 


20 


1 


500 


1.28121 


0.001 


11 


1 


500 


1.28107 


0.001 


25 


5 


500 


1.28107 


0.001 


20 


1 


1 


1.285xx 



TlR 


T"UV 


7"mid 


^■imposed 


2 

m 


0.001 


20 


1 


500 


4.12131 


0.2 


20 


1 


500 


4.12129 


0.001 


11 


1 


500 


4.12131 


0.001 


25 


5 


500 


4.12131 


0.001 


20 


1 


1 


4.12131 


0.001 


20 


1 


500 


6.59516 


0.2 


20 


1 


500 


6.59519 


0.001 


11 


1 


500 


6.59516 


0.001 


25 


5 


500 


6.59516 


0.001 


20 


1 


1 


6.5951x 



Table 6: A few sample mass values m 2 and their variation with some of the program 
parameters. ni mpose d is the number of non-adaptive, imposed steps, within which 
we allow for adaptivity. "xx" means the plot lacks detail to resolve further digits. 
The conclusion is that even under strong changes of program parameters, roughly 
three significant digits remain stable. 



with double precision, e.g. to generate the plots in fig. HI the C code takes about one hour 
on a 3 GHz Pentium 4 PC running Linux. 

Many of the fields vary quickly in this system, and it turns out that the stepsize adjust- 
ment in typical adaptive routines does not allow sufficiently fast change, so we imposed a 
certain number of steps ni mpose d that are not adaptive, and then allowed adaptivity within 
those. Those steps were exponentially spaced, 



/ TtjV \ "imposed , 

Ti+1 — I I T i ) 2 — 1, ... , ^imposed • K^- 1 -) 

\ t irJ 

So ^-imposed is a program parameter that when set to 1 reduces the numerical integration 
to standard Runge-Kutta 4th-5th order over the entire range [ttr, ttjv]. 

We performed some more extensive tests, but a few sample tests are shown in Table 
O along with the variation in the answer, the mass value m 2 . 

So far, this all referred to the midpoint determinant method. We also tried some other 
approaches, as outlined in Table [71 Each of the other approaches has some redeeming 
feature, but for computing the KS spectrum, the midpoint method was superior. As can 
be seen from the table, the main challenge in using the midpoint method is that one 
needs to find approximate asymptotic analytical solution in both the UV and the IR, but 
this can be done using symbolic computation^ We would expect that for backgrounds 
more complicated than KS, the midpoint method will remain the most efficient of those 
in Table [7J 
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For example, Maple can output optimized Fortran or C code implementing those solutions. 
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Method 


Stable 


Multifield 


Avoids Exact 


No fitting 


Standard shooting 


X 


X 


UV 


/ 


Blind shooting [47, App.B] 


X 


X 


/ 


/ 


UV — > IR Determinant 


X 


/ 


IR 


X 


IR — > UV Determinant 


x// 


/ 


UV 


X 


Midpoint Determinant 


/ 


/ 


X 


/ 



Table 7: Sketch of numerical approaches. Stable means numerical stability, in the 
empirical sense of Table El for the examples we studied. Multifield means that it 
generalizes easily to systems with many fields. Avoids Exact means one does not 
need exact analytical solution of approximations of the equations of motion in the 
UV or IR region, or both. No fitting means one does not have to fit the numerical 
solutions to prescribed analytical behavior. 
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